On the quantum (in) stability in cavity QED 
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The stability and instability of quantum motion is studied in the context of cavity quantum 
electrodynamics (QED). It is shown that the Jaynes-Cummings dynamics can be unstable in the 
regime of chaotic walking of an atom in the quantized field of a standing wave in the absence of any 
other interaction with environment. This quantum instability manifests itself in strong variations 
of quantum purity and entropy and in exponential sensitivity of fidelity of quantum states to small 
variations in the atom-field detuning. It is quantified in terms of the respective classical maximal 
Lyapunov exponent that can be estimated in appropriate in-out experiments. 
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The problem of stability of quantum dynamics has at- 
tracted a great interest by its own right and in relation to 
the field of quantum information and computation. Clas- 
sical instability is usually defined as an exponential sep- 
aration of two nearly trajectories in time with an asymp- 
totic rate given by the maximal Lyapunov exponent A. 
Perfectly isolated quantum systems are unitary and can- 
not be unstable in this sense even if their classical limits 
are chaotic pj. It is well known that quantum coher- 
ence is destroyed due to interaction with an environment 
0, Q which is usually modeled by a heat bath with in- 
finitely many degrees of freedom. Environment-induced 
decoherence causes quantum-entropy increase which is 
dominated by the classical Lyapunov exponents |4(. In 
an alternative approach 0, |(| the quantum instability is 
proposed to be measured by the decay of the fidelity or 
the overlap 



/(r) = |<*i(r)|* a (r))| 



(1) 



of two wave functions and ^2, identical at r = 0, that 
evolve under slightly different Hamiltonians. 

In a number of numerical studies 0, H, 0, E3 for a va- 
riety of classically chaotic models it has been established 
that the overlap decay may be algebraic, Gaussian, and 
exponential. The strength of perturbations in Hamilto- 
nians and other factors determine which of these regimes 
prevails. 

In this letter we show that instability of quantum 
dynamics and its exponential sensitivity to initial con- 
ditions and small variations in parameters may occur 
in a paradigmic cavity-QED system with a single envi- 
ronmental degree of freedom. To specify the problem 
we consider the standard model in cavity QED Jaynes- 
Cummings Hamiltonian pd| 
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— hQo (a* a- + a<7+) cos kf£, (2) 

which describes the interaction between a two-level atom 
with lower, |1), and upper, |2), states, the transition fre- 
quency uj ai and the Pauli operators a± z and a quantized 



electromagnetic-field mode with creation, a', and anni- 
hilation, a operators forming a standing wave with the 
frequency uif and the wave vector kf in an ideal cavity. 
The atom and field become dynamically entangled by 
their interaction with the state of the combined system 
after the interaction time t 



|*(t)) = £> n (i)|2,n) + 6„(t)|l, 



(3) 



to be expanded over the Fock field states \n), n — 
0,1,.... Here a n (t) = a n (t) + i(3 n (t) and b n (t) = 
Pn (t) + irj n (t) are the complex- valued probability ampli- 
tudes to find the field in the state |n) and the atom in 
the states |2) and |1), respectively. In the process of 
emitting and absorbing photons, atoms not only change 
their internal electronic states but their external trans- 
lational states change as well due to the photon recoil 
effect. If atoms are not too cold and their average mo- 
menta are large as compared to the photon momentum 
hk f , one can describe the translational degree of freedom 
classically. The whole dynamics is now governed by the 
Hamilton-Schrodinger equations 01 that have the fol- 
lowing normalized form in the frame rotating with the 
frequency u>f(n+ 1/2): 

X = Kp, 

00 

p = -2 sinx^ Vn + 1 (a n p n+1 + /3 n rj n+1 ) , 

n=0 



--j3 n - \Jn + lrfo+i cos a:, 



(4) 



$n = „Ot n + y/n + lPn+1 COSX, 



Pn+1 
Vn+1 

where x 



-7? n+ l - \/n+ IfinCOSX, 



- 2 Pn+1 + v 7 "- + la„ cos x, 



kf (x) andp = (p) /hkf are the atomic center- 
of-mass position and momentum, respectively. Dot de- 
notes differentiation with respect to dimensionless time 



2 



t = Clot, where Oq is the amplitude coupling constant. 
The normalized recoil frequency, k = hk^/m a flo -C 1, 
and the atom-held detuning, S = (ujf — uj a )/Clo, are the 
control parameters. 

Inspite of existence of an infinite number of the inte- 
grals of motion 



R n = al + Pi + Pl+i + vl+i = const, Rn ^ 1 ( 5 ) 

n=0 

and conservation of the total energy 



n=0 

CO 

-2cosa:y^\/n+l (a K p Tt+ i + 0nVn+i) , (6) 

n=0 

the Hamilton-Schrodinger system 10} is, in general, non- 
integrable. The type of the center-of-mass motion de- 
pends strongly on the values of the detuning 5. In the 
limit of zero detuning and with initially excited or deex- 
cited atoms, the optical potential disappears, and atoms 
move with a constant velocity x = npo. The quantum 
evolution is periodic with the period ir/npo, and exact 
solutions for purity, von Neumann entropy, fidelity /(r), 
and other quantum characteristics can be found in the 
explicit forms. For example, the atomic population in- 
version at 6 — is the following: 



z(r) 
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z n — a n~l~ Pn~ Pn+1 ~ Vn+l' 



With the detuning being large, \5\ 3> 0, the optical po- 
tential is shallow, atom moves with almost a constant ve- 
locity, ~ npo, slightly modulated by the standing wave, 
and its inversion oscillates with a small depth (except- 
ing for the case of the so-called Doppler-Rabi resonance 
with maximal Rabi oscillations that occur at the condi- 
tion | <5 1 = Kpo [l3|). If the atomic kinetic energy, np 2 /2, 
is not enough to overcome barriers of the optical poten- 
tial, the atomic center of mass oscillates in one of the 
potential wells. 

Numerical simulation shows that there exist condi- 
tions, defined mainly by the values of the detuning, when 
atoms move chaotically in a cavity. This type of motion 
may be called a chaotic or random walking, and it is 
quantified by positive values of the maximal Lyapunov 
exponent A. In Fig. \I\we show by the dotted line the de- 
pendence A(i5) computed with Eqs. (QJ and the following 
initial conditions: xq — 0, po — 25, the atom is prepared 
in the state |2) and the field is initially in a coherent 
state with the average number of photons n = 10. The 
normalized recoil frequency is chosen to be k — 0.001, a 
reasonable value with usual atoms in a high-quality op- 
tical microcavity in the strong-coupling limit. Stability 




FIG. 1: Quantum-classical correlation between the dependen- 
cies of the variance of quantum purity, op and the maximal 
Lyapunov exponent A (in units of CIq) on the atom-field de- 
tuning 8 (in units of S7o). 



of the computation with respect to truncating the set (@J 
was checked. In most the cases n = 100 was taken. 

The entanglement between the internal atomic and 
field degrees of freedom can be characterized by the quan- 
tity known as purity 



P(r)=Tr ^ (r ) 5 
where p a {T) is the reduced atomic density matrix 



Pa{r) 



(n\p(T) \n) 



(8) 



(9) 



71 = 



with the total density matrix to be p(r) = |^(r)) (^(t)\. 
Purity is maximal if an atom is in one of its energetic 
states |1) or |2); i. e. P max = Tr a p\ = Tr a p a = 1. Purity 
is minimal if p a = 1/2, i. e. P m i n = 1/2, where / is the 
identity matrix. In terms of the probability amplitudes, 
it is given by 



\n=0 / \n=0 ) 

/ OO \ 2 / OO 

"2 ^ ( a nPn + finVn) ~2 £ (a n T] n + (3 n p„ 



\n=0 



vn=0 



(10) 



The root mean square variance of purity, op = 



(P 2 ) — (P) has been computed in the range of the 
detuning |<5| ^ 2 at the same conditions as it was done in 
computing the maximal Lyapunov exponent A. Irregular 
oscillations of ap occurs on the same interval of | <5 1 < 1, 
where A > (see Fig. Computing the von Neumann 
entropy, S = — Tr a (p a lnp a ), we have found the same 
correlations of its variance with A. 

The chaotic centre-of-mass walking has fractal proper- 
ties. Placing atoms at the point x = with the same 
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FIG. 2: (a) Fractal set of the initial momenta po (in units 
of hkf) of atoms that leave a one- wave length cavity after 
m turns, (b) Sensitive dependence of the atomic position 
x (in units kj 1 ) on the initial momentum po- (c) Sensitive 
dependence of the output values of the atomic population 
inversion z out on its initial values Z[ n . Control parameters 
S = 0.4 and k = 0.001. 



initial conditions and parameters but with different val- 
ues of initial momenta po , we compute the time T(po), 
the atom with a given value po needs to reach one of the 
nodes of the standing wave at x = —ir/2 and x = 37r/2, 
and the number of times, to, when it changes its direction 
of motion. The scattering function T(po) is found to have 
a self-similar structure with singularities on a Cantor-like 
set of initial values of momenta pq. In Fig.|2K we demon- 
strate the mechanism of generating this set at 5 = 0.4. 
There are two sets of atomic trajectories with T — > oo, 
the countable one consisting of separatrix-like trajecto- 
ries, corresponding to the ends of the intervals in Fig. Et-i 
and the uncountable one consisting of trajectories with 
to = oo. The chaotic motion can be, in principle, verified 



in experiments on 1 D-scattering of atoms at the stand- 
ing wave. Fig. |2t> demonstrates sensitive dependence of 
the atomic positions on po at a fixed time moment. A 
smooth segment of this function in the range po < 20 
should be attributed to atomic oscillations in the first 
well of the optical potential since these values of po are 
not enough to overcome the respective potential barrier. 
When po exceeds a critical values, atoms leave the well, 
and it is practically impossible to predict even the sign 
of the atomic position. The so-called predictability hori- 
zon can be estimated as follows: t p ~ A^ 1 In (Ax / Axq) , 
where Ax is the confidence interval and Axq the initial 
inaccuracy in preparing initial atomic positions. In order 
to demonstrate the quantum-classical correspondence in 
the chaotic regime we compute the dependence of the 
output values of the atomic population inversion z out at 
a fixed moment on its initial values z; n with the other 
conditions and parameters to be the same. Fig. illus- 
trates that this function is regular in the range |zi n | < 0.5, 
where the atomic center-of-mass motion is regular, and is 
irregular if 0.5 < |zi„| ^ 1 where atoms move chaotically. 

To quantify instability of quantum evolution in cavity 
QED we compute the decay of the fidelity /(r) which is 
the overlap Q of two states |^i(r)) and (t)), iden- 
tical at r = 0, that evolve under two Hamiltonians @ 
with the slightly different detunings A<5 = Si — 6. 
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FIG. 3: Decay of the fidelity of quantum motion 1 — f(r) 
(logarithmic scale) in the chaotic (thick and thin lines) and 
regular (dotted line) regimes. Time r is in units of . 

We have found previously (see Fig. 0) that with the 
initial momentum po = 25 and 8 = 0.4 the type of 
atomic motion depends strongly on the initial atomic in- 
version z(0). If an atom is prepared initially in one of 
its energetic states, i. e. z(0) = ±1, its classical and 
quantum dynamics are unstable, whereas they are sta- 
ble with z(0) = under the same other conditions. In 
Fig. El we show for convenience the decay of the quantity 
log 10 (l — /) in the regimes of chaotic walking (thick and 
thin lines, z(0) = ±1) and regular motion (dotted line, 
z(0) = 0) with AS — 10~ 4 . In the chaotic regime the 
fidelity decays exponentially with the rate A ~ 0.04 to 
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be equal to the maximal Lyapunov exponent computed 
with the set @. This result does not depend on the 
values of differences in the control parameters AS. The 
fidelity practically does not decay in the regular regime 
at z(0) = 0, and the respective maximal Lyapunov expo- 
nent was computed to be zero. 

We emphasize that sensitive dependence of quantum 
motion both on initial states and parameters may arise 
with an atom in a quantized cavity field. Single chaotic 
degree of freedom, arising naturally when we take into ac- 
count photon recoils, provides quantum instability and 
irreversibility. We do not need an infinite bath or any 
kind of noise for that. The quantum instability has 



been shown to be quantified by the respective maximal 
Lyapunov exponent providing a quantum-classical corre- 
spondence. 
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